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Abstract. In this work we study, under the Stratonovich definition, the problem 
of the damped oscillatory massive particle subject to a heterogeneous Poisson noise 
characterised by a rate of events, A(i), and a magnitude, 3>, following an exponential 
distribution. We tackle the problem by performing exact time-averages over the noise 
in a similar way to previous works analysing the problem of the Brownian particle. 
From this procedure we obtain the long-term equilibrium distributions of position and 
velocity as well as analytical asymptotic expressions for the injection and dissipation 
of energy terms. Considerations on the emergence of stochastic resonance in this type 
of system are also set forth. 
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1. Introduction 

Stochastic processes liave long surpassed ttie limits of a mere probability theory subject 
establishing itself as a topic of capital importance in various disciplines which go from 
molecular motors to internet traffic analysis [H E] [3]. In what concerns statistical 
mechanics, we owe its introduction to the study of Brownian motion by means of the 
Langevin equation jl]. Explicitly, a random term was added to the classical laws of 
movement in order to emulate the collisions between the particle under study and the 
particles of the supporting medium assuming that the jolts between particles cannot 
be deterministically written down [Sj E]. This approach was promptly adopted in 
other problems, particularly those coping with systems that exhibit a large number 
of degrees of freedom. Therehence, the allocation to the noise of the microscopic 
features and playing interactions of the system has become a common practice. Due 
to the accurate descriptions of diverse systems, the use of the Wiener process (driftless 
Brownian motion) [3, |8] has assumed a leading role in widespread fields. However, this 
corresponds to a very specific sub-class in the wider Levy class of stochastic process [9]. 
Specifically, the Wiener process corresponds to the case of a stochastic process for which 
the increments are independent, stationary (in the sense that the distribution of any 
increment depends only on the time interval between events) associated with a Gaussian 
distribution thus having all statistical moments finite. Moreover, its intimate relation 
to the Fokker-Planck Equation has put the Wiener process and its adaptive processes 
in the limelight. 

Despite the ubiquity of the Wiener process, several other processes, either in Nature 
or man made, are quite distinct [lOl [TTl [12]. Particularly, they can be associated with 
(compound) Poisson processes in which the independent increments befall upon a rate 
A (t) with its magnitude related to a certain probability density function. For this kind 
of noise, it was previously shown that the traditional Fokker-Planck approach cannot be 
applied because the Kramers-Moyal moments from the third order on do not vanish [13] . 
Consequently, the equation for the evolution of the probability density function cannot 
be exactly written as a second-order differential equation, but it maintains its full 
(infinite) Kramers-Moyal form instead [13]. As it turns out, a complete solution is 
quite demanding for most of the cases. 

In this manuscript we study a damped harmonic oscillator subjected to a random 
force described by a heterogeneous compound Poisson process which can be replicated 
by means of a RLC circuit with random injections of power, or several other dynamical 
processes at the so-called complex system level. Moreover, such kind of problems often 
conduce to the emergence of stochastic resonance [T^ 115] that is found in so diverse 
problems spanning from neuroscience and Parkinson's disease [16] to micromechano- 
electronics [H]. Although this specific system (and its variants) has been studied 
by different authors [1], in this manuscript we survey the problem assuming a very 
fundamental and different approach. Explicitly, we make direct averages over the noise 
for different quantities of interest, in the reciprocal space-time (Fourier-Laplace) domain. 
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by a method that can be apphed to the treatment of coloured noise [TH], multiple 
types of noise [19], thermal conductance problems [20], or work fluctuation theorems for 
small mechanical models [21] . Namely, within the Stratonovich definition, we first focus 
our efforts on the evaluation of the steady state probability density function (PDF). 
Afterwards, we study the evolution of the injected (dissipated) power into (out of) 
the system as well as the total energy of the system and the emergence of stochastic 
resonance. 



2. Exactly solvable model 



Our system can be described by two coupled stochastic differential equations: 

x{t)=v{t), (1) 

Mv{t) = -kox{t) - -fv{t) + rj(t), (2) 
where ?7 is a compound heterogeneous Poisson noise, 

v{t) = J2,'^it)Ht-ti), (3) 

for which the shots (events) occur at a rate A (t). Specifically, we will study the following 
time dependency, 

A(t) = \o[l + A cos{uj t)] , (0 < A < 1) . (4) 

The case A = yields the standard homogeneous case. We have opted for a sinusoidal 
dependence of the rate A (t) because of its ubiquity in a wide range of phenomena [22] . 

The Poisson process, P, is a continuous-time stochastic process belonging to the 
class of independently distributed random variables whose generating function between 
t and t + At is 

"t+At 



Gt,t+At (z) = exp 



A (t') dt' 



which is defined by the rate of events, A (t), such that. 



{nit)) 



A (t') dt', 



(5) 



(6) 



(we use the notation ((...)) to represent averages over samples whereas (...) symbolise 
averages over time). According to its definition, it is only possible to have a single event 
at instant t. Therefore, from Eq. ([5]) and considering the limit At — >■ 0, one can compute 



the probability that an event occurs is given by A (t) dt, while 1 — A (t) dt 

is the probability of zero events. Again, from Eq. ([5]), we can analyse the probability 
distribution of the inter-arrival time PDF, T(r). The probability that having an event 
at time t the next event occurs at t + r is given by. 



p{T\t) 



d_ 

dr 



exp 



t+T 



A {t') dt' 



(7) 
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For a heterogeneous Poisson process, one must still pay attention that each instant has 
a different weight in the calculation, exactly because of the time dependence of A. This 
weight, / (t), is defined by 

f{t)=\{t)/ r \{t')dt\ (8) 



for a heterogeneous Poisson process taking place in the time interval from to T. 
Thence, combining Eqs. ([7]) and (jH]) and integrating over time one obtains, 

r(r)= r ^ f{t)p{T\t)dt. (9) 



Concerning the amplitude, $, despite being possible to consider several distribution 
functions, we will restrict our study to the classical exponential probability density 
function for white shot-noises. 



P ($) = $-1 exp 



whose n-th order raw moment is, = 71.!$". 

Going back to Eq. ([3]), we can define the Poisson patch, /, between t and t + r, 

rt+T rt+T 

I (t) =F{t + T) - F{t) = dr= 7] {t') dt', 

Jt Jt 

from which we can set out, 

dl (r') = J dF = j 7] {t') dt', 

where we have omitted the time dependence of the patch for the sake of simplicity. The 
average of the Poisson patch is, 

a(r)),= (/(r)) = <i'^'^^A(t') dt'. 

and its covariance, 

/t+Tl /•t+T2 
V it') dt' V in dt". 

which yields after averaging, 

Xit") 5it" -t') dt"dt' 

'•t+min(T2,Ti) 



/t+miiHT2,Tij 
Xit') dt'. 



These moments can be straightforwardly generalised to, 

dti...j^ dtnVitl)---Vitn). 

The noise cumulant averages are, 

pt+At pt+At 
{I it)""), = dt,... dtn ivitl) . . . Vitn))c = A it) WAt. 
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This implies that the noise cumulant correlations are identified as 

(r/(ti) . . . 7]{tn))c = KhW^ih - ts) . . . 5{tn-l ' tn) ■ 

Accordingly, 



(10) 



t+At 



l[dU{r]{t^)...v{tn))c 



i=l 



t+At 



n-1 



Yl dti Ao (1 + A cos(a; ti)) JJ 5{tj 

i=i j=i 

ft+At 

Ao$" / dti (1 + Acos{uti)) 



Ao 



A / At , , . 

At + 2— sm u cosfcj t) 

CO \ 2 



Taking into account the limit, At — > 0, the previous expression tends to 
Ao$" At(l + Acos(a;t)) and thus, 

/t+At pt+At 
dti...j^ dtn ivik) ■■ ■vitn))c = >^o'^['^ + A cos (ut)] At = X{t)W At. (12) 

In the case n = 1, Eq. f ll2p satisfies the relation (/ (t)) = X (t) $ At, as expected. 

Throughout this manuscript we employ the Stratonovich representation for the 
noise, 



A (t) 



At 



/~\Jt 2 - 2 

where the effective noise in a interval dt is computed as the average of the noise at t 
and t + dt. 



2.1. Laplace transformations 

Taking the Laplace transformations of Eqs. ([T]) and ^ (with Re(s) > 0) we obtair(||, 
sx{s) = v{s). (13) 

Defining, 

^^'^ ^ + + ^ = " " ^^^^ 

we can write. 



xis] 



MR{sy 

the singularities of which (the zeroes of R{s)) are located at, 
K± = —9 ± i fi. 



(15) 



(16) 



with e = = I? and Q 



\ In the Laplace transform we have assumed that both the initial position and the initial velocity are 
equal to zero. We are interested in asymptotic effects and the initial memory terms vanish in that limit. 
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For the Poisson process with time- dependent rate, Eq. (jlj), the Laplace transform 
of the noise averages yields, 



{fl{Zi) . ..fl{Zn)) 



/•oo " f " 1 

c= / TT citi exp <^ - V'zit^ W?7(ti) .. .?7(t„)), 

-^0 ^=l I 1=1 J 

poo " 

= XoW T\dti 6iti -t2)... 5(t„_l - tn) 

•Jo ,-1 



i=l 

X [l + Acos(wti)]exp|-^2;iti|, (17) 

for which we can separate out its homogeneous and heterogeneous parts. For the former, 
we obtain, 



<>oo r _ 

Xi = Ao / dti exp <^ -ti V Zl 

-^0 I 1=1 



^ Aq 

En 1 
1=1 

whereas the latter is given by. 



/■oo r _^ _ 

l2 = XqWA / dti cos{uti) exp <^ -ti ^ 



= Ao f ^n ^ ■ + \ . ) • (19) 

The Laplace transform for the noise cumulants is thus given by: 
(f/(2;i)...r/(2;„))^ = Xi+X2. (20) 

3. Averaged steady state 

Instead of equilibrium conditions, a periodically forced system reaches a periodically 
driven state that is characterised by periodic variations of the averages and cumulants 
of its variables [23]. This behaviour can now be explicitly obtained by the method 
previously mentioned [18], that we shall use to study the Poisson process. However, 
taking the time average for the distribution function does make sense given that 
important quantities, such as the injected and dissipated energies, are well understood 
when represented by their time averaged values. In the following, we develop the 
techniques needed for the exact solution of Eqs. and ([2]) at long times. Our results are 
comparable with those obtained from the analysis of a single (and sufficiently long) run 
of the process in which all the values of the observable are treated as equally distributed. 

We average the probability distribution over time to obtain the cumulant 
expansion [H [TBI EH] 

Pss{x,v) = V / ^ e^(Q^+P^) V V ^ ^n,,,™ 
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exp 

ZTT ZTT ^-^ n\ ml 



n,m=0:(r?i+n>0) 

Let us work out the exact form of the cumulants in terms of the Laplace transforms 
of the noise. Applying our definition for computing averages we have, 



(x"^;™) = lim^ / dte-'^ 

Jo 



limlim pTT^pn^ 



where the integration path for the variables in complex space is the same as in 
Refs. [18l[19]. 

4. General steady state distribution 

Following the previous section, the Poisson steady state distribution exactly yields. 



(21) 



/ ^ iTT ivr ^-^ n\ ml 

°° \^n+m=0:{n+m>0) 

We may then split Eq. (l2Ti) into two parts: term V^^^ which arises from the time- 
independent contribution, 

i-oo ,y 2vrn™r'/?(ij>.+6) i?(-Eri""'(ip.+e))' ^ ^ 

and the remaining terms arising from the periodic forcing term, 

pW=limM r^^Tl'dP'^ ' U:=,iiPk + e) 1 

""^ 2 y_ n 2vr.±ia;n:;r'^(iP. + ^)^(±i^-E:;r(iP. + ^))' 

Unrolling Eq. ( l22l) . we can determine the different contributions that are allowed to 
emerge by taking into consideration the different powers of x and v in the cumulant. 
Accordingly, term \E^i represents cumulants of zeroth order in the velocity, (x"")^, 

E f " - M ,. jxr: (-) 



j=0 



and in the position, (w™)^ 



' c' 

m— 1 

m 



^ _ V- I - M i"^ (-1)""^' kj k--^-' [j , (m - J - 1) ] 

^' / [l,-l]"^-Mj,(m-j)][(j + l),(^-J-l)]' 



i=o 
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Term '^2 describes cumulants which are such as (x*^ ^v)^ 



n-2 
j=0 



n-2 
J 



[j, {n-j-2)] 



[1,-1]""' [j, (n-j-l)][(j + l) ,(n-j-2)] 



(25) 



and last term \l/3 represents the remaining combinations of powers of and f ™ with 
I + m = n and / > 2, 

m n—m—l 



^3 M 



E E 

i=o «=o 



n — m — l 
I 

1 



(-1)"-'"+^ [^■^(^_^-)] 



[1,-1 



in— m— 1 



[1,-ir 



X 



[U + I + 1) ,{n - m - J - I - + I) ,{n - m - J - l)y 

in which we have the used the following curtailed notation, 

(a + 6 /€_) = [a, b] . 

Therefore, bearing in mind the last definition of the cumulant part of the probability 
distribution, Eq. f l22|) . we can write it as: 



exp ^ Ao(n + m)! 

n,m=0:(m+n>0) 



qn pm 

n\ ml 

n+m—l 



(i$) 



X 



m=2 



(27) 



where ^-"(3;,^) [/ (Q, P)] represents the 2-dimensional Fourier Transform into position- 
velocity real space, (x,w). 

Owing to the factor z/ {z ^lu) in terms p(^) we are able to verify that after 
performing integrations following the appropriate contour, the z terms hold on up to 
the end, so that when we finally compute the limit of 2; — ?■ both contributions vanish. 
Therefore, for this type of heterogeneity, the steady state distribution bears out the 
same result as the homogeneous Poisson process with constant rate Aq. This is quite 
understandable since we are making a long time average in which the contribution of the 
periods whose rate is larger than Aq kills off the contribution arising from the periods 
in which the rate is smaller than Aq because of the symmetry of the rate around Aq. 

Regarding the marginal steady state distributions. 



Pss {x) = / Pssix,v) dv, 



and 



Pss (v) = / Pss{x,v)dx, 



(28) 



(29) 



we start with the probability distribution of the position and following our procedure 
we obtain, 



Pss (•^) J'x 



exp 



la; 



. n>0 



(30) 



X 



(26) 
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whence we can identify the cumulants, 

«:„=^ = n!Ao (i<l)"*i.. (31) 



Using the property of Pascal triangles, 

n — 1 \ I n — 1 




(32) 



we can write the sum in Eq. fl23|) over the index j as, 

n— 1 

5^ [...] = n! [1, -1]"-^ — , (33) 



where. 



Dn = \{[3,n-3]. (34) 
i=o 

Accordingly, the cumulants of Pss{x) are 
Allowing for Eq. (135!) . we explicit the average. 



fco 

and the second-order moment, 

■2 -o A 



2 ^0 



We have implemented a computational procedure to numerically compute the 
probability density function of the position at the steady state. Our exhibited numerical 
results were obtained for different values of M, ko and 7 and fixed values of Aq = 10, 
$ = 1 and u = TT following the implementation described in Appendix A From cases 
A and B (see values in the figure caption), we can understand that the mass does not 
impact in both the average and standard deviation and that B and C tally, although each 
is based on a homogeneous and a heterogeneous process, respectively. Comparing cases 
A, B and C we verify that the lighter the particle, the skewer the distribution {x). 
The results of cases B, C and E help us show that 7 does not affect the average and 
finally case F sketches the influence of k^. Case D permit us to follow the dependence on 
the mass. Lighter particles have a skewer distribution. The results are clearly different 
if we consider a symmetric noise with $ = 0. As we will see in the next section, the 
positive total injection of power associated with this case is replaced by a zero average 
injection of power related to the steady state average position being zero. 

The picture is very much the same for the marginal distribution of the velocity. 
Namely, from Eq. (I29l) we get. 



Pss {v) = J^v 



exp <^ ^ Ao P"* 



. m>0 
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-5 5 10 15 20 25 30 35 



Figure 1. Numerically obtained probability density function pss {x) vs position x for 
various cases with Aq = 10, $ = 1 and the noise defined by Eq. ^ with u = tt. 
Following the legend in the figure we have the respective cases, A: M = 1, fco = 1, 7 = 
I, A = 0, B: M = 10,fco = 1,7 = 1,A = 0, C: M = 10,fco = 1,7 = I, A = 1/2, 
D: M = 0.1, fco = 1,7 = 1,A = 0, E: A/ = l,feo = 1,7 = 2,A = and F: 
M= l,fco = 10,7= 1,A = 0. 



and consequently the cumulants are, 

= (^ = m!Ao {i^r^i,. (36) 

For the cumulants of the marginal velocity distribution, the calculation turns out 
much harder and haplessly we have not managed to write it in a compact form as 
Eq. (|30|1 . Nevertheless, we can still write some of them explicitly, such as the first, 

and the second cumulants, 

^ ^ ' M7 

In Fig. m we plot the results of pss {v) for the same the numerical implementation of 
Fig. [H Once more, we can understand the independence of the probability distribution 
regarding the amplitude A in Eq. (j4]). We can also notice the influence of the mass, 
lighter particles are more sensitive to the noise and thus, for the same noise intensity, 
they achieve larger positive values of the velocity. Moreover, it is visible that pss {v) 
might be strongly positively skew for light particles. On the other hand, if we consider 
an average over samples the effect of the amplitude A and frequency u of the Poisson 
noise. 



5. Injection and dissipation of Energy 




Figure 2. Numerically obtained probability density function pss {v) vs scalar velocity 
V for the same parameters sets of Fig. [T] 

An interesting element of study, especially in practical applications such as presented 
in Ref. [T71 [T3], concerns the time evolution of energy (related) quantities. This can be 
checked by heeding the fact that variations of energy in an isolated system equals the 
total work done by the external forces acting on it, in this case the fluctuating force rj 
and the dissipative one — 7f : 

J2 Wf.... = (37) 

f 1 2 1 2 

rjdx — -y / vdx = -M v (t) + -ko x (t) 

The injection of the energy balance can be also analysed by determining the 
evolution of the total energy of the particle, 

^ [V (t) V (t) - [tf] dt = \m V {tf\'Zl + \ko X {t)XJ, ■ (38) 

In what follows, we shall omit transient terms, keeping the more interesting 
asymptotic ones whenever possible. 

5.1. Energetic considerations 

The average values of [^'(t)]^ and [a;(t)]^ that we use in the computation of the total 
energy can be obtained once more using the Laplace representation. 



V {tf = I I V (ti) V (ta) 5 (t - ti) S{t- ta) dti dt 




2 



liin [ [ e(i«i+''?^+2.)* ^ (i + ^) ^ (i + ^) ^ 

e->0 / / ZTT ZTT 



whence by averaging and taking into account the second cumulant definition, 
lim — - 
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{fi (i qi + e) fi (i q2 + e))^ dqi dq2 
R{iqi+e) R{iq2 + e) 2ti 2n' 

and, 

{^x[t)) {x[t)) -lim^^yye R{iq,+e)R{iq, + e) 2n 

where the asymptotic solutions are, 

/ _/ m\2 _>^^oM^i^'d'+d^o')cos{cut)+co^9^sm{ujt)]Xo^' , 



^ A [2 ^ (w^ - 3 wo^c^^) COS (w t) - 3 wo^w^ sin (w t)] Aq^^ 

M2 (4 ^2 + cj2) (^(^2 _ ^2)2 ^ 4 ^2^2 j 

A (cj^ + 8 wo'^w ) Ao$^ sin (w t) 



2 



M2 (4 ^2 + c^2) (^(^2 _ ^2)2 ^ 4 ^2^2 



and 



(r(f?\ -(r(f))' f 2 Ao 4a; (a;2 _ 4 a;o2 - 8 ^2) A l>2 Ap sin (a; t) 

8 (3 ^ U;2 _ 2 ^ t^g2) ^ ^2^^ ^) 



4 



1 Uo'^U^ + 16 ^2cj2 + 16 cOo^) (4 ^2 + ^2) JVf 2 ' 



We must now take into account the square values of {v (t))^^^ and {x {t))^^^, which, 
for all times, are given by, 

^ ^ e^oM J 271 Riiqi + e) 



e-i>o M 7 27ri?(igi+e)l^igi + e 2\igi + e — io; igi + e + icu 



and, 



e^o MJ 27r_R(igi+e)[ 2\igi + e — iw igi + e+iw 
After integrating over the poles above, we obtain the asymptotic behaviour 

A [(u;2 _ ^q2-j g-j^ (a; t) + 2 w ^ cos (w t)] u Ao$ 



(x(t)) 



cc;2-a;2)- + 4^2 ^2 
$ Ao v4 [(a;2 - wq^) cos (w t) - 2 61 a; sin {u t)] Ao$ 



^2 jyj 



[w2_^2)2^4^2 ^2 



where it can be easily seen that. 

We observe that the results above are those expected when we interpret the 
oscillating rate Poisson process as a periodic forcing acting upon a damped harmonic 
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oscillator. As expected, the amplitude of motion shows the typical resonant behaviour. 
In order to illustrate this behaviour, we plot in Fig. [3] the quantity. 



2 

asy 



\ 



(39) 



as a function of the amplitude of the oscillating contribution. A, and the frequency of 
these oscillations, u. In the both planels is evident the emergence of a maximum at a 
frequency of the heterogeneous Poisson rate A {t) equal to ^Juj^ — 29"^. 




^01234 

w 



Figure 3. Upper panel: y{x{t))^gy vs amplitude of the heterogeneous part, A, 
and frequency, w, with M ~ 1, "f ~ 1, ko ~ 1, Xq ^ 10, $ = 1. The maximum, 
characterising a stochastic resonance phenomenon, occurs at uj^cs = 

ytl^^~~202 that 

defines the grey dashed line in the lower panel containing the plane cuts for A — 
(black line), A — ^ (dashed red line) and A = 1 (dotted green line). For this case the 
maximum occurs at 1/V2. 

Adding all the terms and taking the limit t — > oo, we obtain the equilibrium energy 
of the system. The energy is composed of an oscillating term, with time zero average, 
and a constant term El,j: 
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$2 

7 



2Mul 



4M 



14 

(40) 



It is worth mentioning that in this case we have made exphcit the asymptotic time 
dependence so that our averages are computed over samples and not over time in a 
single sample as we have made in the previous section. 



5.2. Power considerations 

Going back to Eq. (1551) . we can define the two following quantities, 

Ji = v{t)rj{t), (41) 

and 

JD = ~lv\t). (42) 

Physically, both rates, J/ and Jd, constitute changes of energy due to the interactions 
with the thermal bath. Within this context, it is particularly important the study of 
the cumulative changes of energy in the system up to a time t = t, namely the injected 
total, 



JiT= I dtv{t)ri{t), 
'o 



and the dissipated total, 
Jdt 



(43) 
(44) 



-7 / dt v\t). 
Jo 

in which we will apply the same Laplace transform operation in order to better handle 
the noise averages [TH [191 [201 [21] . The dissipation of energy fiux can be written as 



J, 



DT 



nT noo POO 

-7 / dt dti5{t-ti) / dt2 5{t-t2)v{ti)r]{t2) 
Jo Jo Jo 

dqi dq2 e('9i+i<?2+2.)r _ ^ 



lim 

e^O 



lim 

e^O 



oo 27r 

oo 27r 
igi + e) 



oo 27r (i gi + i + 2e) 
^ dq2 e(' 31 +'92+2^)^ - 1 
^ 2tt (i gi + i + 2e) 
(i q2 + e 



v{iqi + e)v{iq2 + e), 



MRiiqi + e) M Riiq2 + e 
which becomes, after taking the thermal average 



?7(igi + e)?7(ig2 + e) 



{Jdt) 



27A0 
NP 
1 



lim 



dq 



dq2 e 



(igi+i g2+2e)T 



1 (igi + e) (ig2 + e) 



00 27r y_oo 27r (igi + ig2 + 2e) i?(igi + e) i?(ig2 + e) 
A f I 1 



X 



igi + ig2 + 2e 2\igi + ig2 + 2e-icj igi+ig2 + 2e + ia; 

2 l^igi + e 2\igi + e — iw igi + e 
I A f I 1 

ig2 + 2e 2 \ig2 + e — icu ig2 + e + ia; 

JDToij) + JDMscij). 



X 
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Under the condition t9 ^ 1 (transient terms are negligible), we can explicitly write 
the results for the integration above as a sum of three contributions: a term proportional 
to r, 



M 



an oscillating term, 



D,osc 



j) 



A (3 ^ - 4 ^ o;^ - 2 ^ a;^ - 4 ^3 ^2^2^ ^) 

M^u (4 + (^^4 + ^2^2 ^ ^4 _ g ^2^2) 
A (8 w + 4 ^2a;3 + ^5 _ g ^3^^2) ^2^2^ ^) ^ 

M^U (^2 + 4 e^) + 16 ^2^2 + Ig ^4 _ 8 ^2^2) - 

sin (2 a; r) + ( 



+ 



a;2)cos(2u;r)] <|2a2^2^2 



M [a;2 (^2 _ 2 u;o2 ^ 4 ^2) + ujlf 

^2 [^2 (^4 ^ _2 ^^2) ^ ^^4] ^2^^2^ ^ ^j^^ W t) 



2M [a;2(^2_2c^2^4^2) 



(45) 



and a constant term, 
$2 Ao 



J, 



$2A2 



7 



+ 



AI.2A2 (a;2 - a;2) 



2Ma;2 [a;2 (a;2 _ 2 a;2 ^ 4 ^2) + ^2 



+ 



^2 [[2 a;2 (,^2 + 2 g2) - ^2 ^ g o;^ {6^ - uq")] $2 A2 
2 M [a;2 (^2 _ 2 u;2 ^ 4 ^2) + u;2]2 

where Jdt = Jdto{t) + 

The injection of energy can be written in a similar way. 



(46) 



J, 



IT 



l*T 1*00 1*00 

dt dti6{t-ti) / dt2 6{t-t2)v{ti)7]{t2), 

Jo Jo Jo 



e^O 



dq 



lim / / 



e-i-O 



27r 
dq 



dq2 e('3i+''?2+2^)^ - 1 



lim / / 



2ir 



2tt (i gi + i g2 + 2e) 

dq2 e('9i+i92+2e)r _ ^ 



{;(igi + e)r]{iq2 + e) 



27r (igi + ig2 + 2e) 



M R{i qi 



-Tjiiqi + e)?7(ig2 



where after taking the thermal average gives. 



{Jit) 



M 



e-5>0 



dq 



lim / / 



+ 



27r 
A 



dq2 e 



(i gi+ig2+2e)r 



27r (igi + ig2 + 2e) i?(igi + e) 



1 



+ - 



1 



igi + ig2 + 2e 2 \igi + ig2 + 2e - icu igi 



+ 



X 



Ao 



ig2 



1 A 

i gi + e 2 

1 A 
21 ^ 2" 



: — + : — 

igi + e — iw igi + e + iw 

1 1 



i g2 + 2 e + i a; 

X 



1^2 



i g2 + e + i 



= JiToir) + JiT,oscir) + J, 



ITc- 
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An important part of the injection of energy flux has to be carefully obtained since, 



[JiToir)), 



M 



r lim 

£->0 



dqi 
2tx 

^ dqi 
^ 27r 



dq^ g{igi+ig2+2e)r _ ^ (igi + e) 



, 2n (igi + ig2 + 2e)2 i?(igi + e) 

(igi + e) 

(igi + e - K+)(igi + £-«;„ 



where the last term in the r.h.s. contains the integration over the the upper arch, 
because the reduction lemma is not valid in this case, and using the relations between 
raw moments of $ as well. 

Then finally, we write the contributions for the injection of energy as 



JiToir) = 
and 

JiT,osc{t) 

and 

JiTc 



M 



M \uj 



2^W2 _2a;2 + 4^2) +^4 



Ao$M sin {(jj t) a [A (wq 



Mu 



+ 



cos (2 a; r) + 4 ujq^ cos {u t)] Aq 



4M [^2(^2 _2a;2 



4^2) +^4] 



+ 



A[2oje (A sin (2 a; r) +4 sin (wr)) - A u"^ cos {u t)]\1^ 



^2 



4M [w2 (^2 _ 2^2 + 4^2) 



(47) 



Ma;2 



+ 



M [w2 (^2 



2a;2 



+ 



4 / 2 



12 ^2) +^2 



4^2) +^4^ 

r 2 2 , 
— OOOqUJ + UJ 



4 9' 



4 M b2 ('a;2 



2a;2 + 4^2) 



{41 



For the total energy flux Je = Jit + Jdt it can be easily seen that the linear term 
on r cancel out, while the constant term becomes exactly the average energy of Eq. pUj) . 
The non-oscillating part of the energy reads, 

i>2Ao (2 Mwo^ + Ao7) (w^ + ujo^) A'^^'^Xq^ 



jno 
'J p. 



AM [a;2 (a;2 - 2^2 ^4^2) 



41 ' 



(49) 



coinciding perfectly with Eq. (140|) . The oscillating term Jit,osc+Jdt,osc do not contribute 
if the time average is taken. Thus, the average energy will in the form of an oscillation 
around the mean. 

As expected, in the long-term the magnitudes of the injected and dissipated 
energy fluxes are exactly the same signalling the emergence of equilibrium. All of our 
calculations are compatible with the plots in Fig. H] whereby we depict the evolution of 
the injected and dissipated power total (mechanical) energy of an oscillator following 
our dynamical equations. 
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600 




600 



Figure 4. Upper panel: Total injected and (symmetric) dissipated power, Jit {~Jdt) 
vs time, r, according to the definitions in the legend. Lower panel: Total energy, Em-, 
vs time, r. The dashed (green) line represents the assymptotic limit given by Eq. (|40|) . 
In both cases we have used the following values: il/ = 10, fco = 1, 7 = 1 $ = 1, Ao = 10 
and A ~ 0. n.h.: The difference between the values of Jit and —Jdt in the upper 
panel is exactly equal to the total energy, Em, that is shown in the lower panel and 
equals the theoretical value given by Eq. (j40|) as well. 



6. Concluding remarks 

In this work we have revisited the problem of the damped harmonic oscillator subject to 
a heterogeneous Poisson process. Our approach, which is carried out by averaging over 
the noise in the Fourier-Laplace space, allowed us to obtain the long-term distributions 
of the position and distributions (joint and marginal). Moreover, we have surveyed 
the interaction between the system and the thermal bath by computing the rates of 
energy that are dissipated from the system and injected into it. As expected, after a 
transient time, both rates balance so that the system achieves a steady state. The 
application of time averages over the position and velocity of the massive partible 



On exact time- averages of a massive Poisson particle 18 

300n • ^ • ^ • 1 

J ■ 

200- 

100- ..• 

0-s 



50 100 150 

T 

350 
340 



330- 
320- 
310- 

300 n • ^ • ^ ' ^ • ^ • , 

200 210 220 230 240 250 

T 

T 

Figure 5. Upper panel: Total injected, Jit, and dissipated power, Jdt, vs time, r, 
according to the definitions in the legend for the heterogeneous case 7\f = 10, fco = 1, 
7 = 1 $ = 1, Ao = 10 and A ~ 1/2. Lower panel: Comparison of Jit for this 
heterogeneous process with a homogeneous process {A = 0) under the same remaining 
parameters. The impact of the oscilations in the noise are clearly visible. 

has allowed obtaining the long-term distributions of the two quantities, which is 
independent of the heteregeneous character of the noise. This last feature will only have 
impact when averages over samples, instead of averages of the time, are implemented. 
Notwithstanding, we having been able to find the effect of the heterogeneity of the 
rate of events be the emergence of resonance effects linking the "natural" frequency 
of oscillation and and the frequency of the time-dependent part of the rate of events. 
This impact is also visible when the total injected and dissipated powers have been 
surveyed. Considerations regarding Jarzynski's equalities as well as modifications on 



-lUU- 



A = 





A = 


1/2 
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the inter-event rule of the noise are addressed to future work. 
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Appendix A. Numerical calculations 

In order to solve in a numerical way our main Eqs. ([I]) and ([2]) we have considered a 
trapezoidal approximation which is reminiscent of the Stratanovich approach to noise 
phenomena, 



for small enough dt. 

In respect of the stochastic part [23] , its calculation can be at least made threefold. 
The first one concerns the inter-event time which is given by Eq. ([9]). Accordingly, 
starting from to we would randomly select a certain time interval, following PDF in 
Eq. ([9]) and we would let deterministic equations evolve up to to + when we would 
add the value of the kick, M~^$ (to + 5t) , to At^'^'^* with $ chosen from the specific 
distribution P ($). Despite its accurateness this procedure is not the most hard-headed 
when it comes to simulating heterogeneous Poisson processes, since we are obliged to 
constantly update the distribution in a rather grinding way. 

The second method is a shrewd procedure of carrying out the numerical simulation 
without having to heed to a very tiny value of the mesh 5t or also to problems with a 
very high rate of occurrence of events A (t). In this case we can determine the (expected) 
number of events that take place within 5t, 



and consider that the overall effect of the noise in that time interval is equivalent to 
the occurrence of a single kick the intensity of which is given by the convolution of 




and the same for the deterministic part of the velocity. 




(ti - ta) 



7[i> {t2)+v{ti)] + ko [X (t2)+x(ti)] 

2M 
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n (t) distributions P ($). Bearing in mind that the events are uncorrelated the resulting 
distribution is given by, 

p„ ($) = j-^ 1 [7, [p mi ■ 

However, it must be stressed that this procedure is half-averaged since it already assumes 
the mean number of events in its implementation, thus leaving all the randomness to 
the resulting amplitude of the added noise. Although we have not tested the following 
assertation, we believe that its application reduces the number of samples needed to 
obtain the same dispersion in the sample set. 

The third way corresponds to our main option, particularly for the figures in the 
text. Specifically, we were intentionally careless about optimising computational time, 
as we have preferred a very conservative approach and a very tight grid. Since we did 
not opt for depicting examples with very high event rates, we went ahead by picking 
a random number uniformly distributed between and 1 and compared it with the 
probability of having an event according with the Poisson distribution with parameter 
jt+st ^ ^^f-^ ^ which is similar to Eq. ( fT2l) . If the random number is the smaller of both 
numbers, then a kick takes place and therefore we need to select the noise intensity 
as described for the first approach. In all the cases we have shown 5t = 10~^. The 
distributions were obtained from a total of 5x 10^ records (10^ samples) made at intervals 
of 10~^ time units. To assure equilibrium we have set apart the first 10^ (100 time units) 
logs of each sample. 
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